mm ~ 
亚 音 时 型 前 缘 形状 对 附 面 层 参数 影响 
于 贤 君 ， 朱 宏伟 ， 刘 宝 杰 
(航空 发 动机 气动 热力 重点 实验 室 , 先进 航空 发 动机 协同 创新 中 心 , 北京 航空 航天 大 学 能 源 与 动力 工程 学 院 , 北京 100191) 
摘要 : ”采用 数值 模拟 方法 研究 了 圆 弧 形 、 椭 贺 形 和 曲率 连续 无 吸力 峰 前 缘 对 叶 型 附 面 层 发 展 的 影响 规律 。 通 过 提取 不 同 工 况 
下 叶 面 附 面 层 参数 并 结合 叶 表 附 面 层 速度 型 的 变化 ， 发 现 前 缘 吸 力 峰 对 叶 型 性 能 的 影响 本 质 在 于 吸力 峰 影响 了 附 面 
的 起 始 发 展 状态 。 前 缘 吸力 峰 扩 压 参数 Dspire 过 强 会 导致 附 面 层 起 始 发 展 状态 恶化 ， 造 成 提前 转换 ， 甚 至 出 现 分 离 泡 ， 
叶 表 附 面 层 迅 速 增 厚 ， 叶 型 损失 增 大 ， 可 用 攻 角 范围 减 小 。 研 究 还 表明 ， 所 设计 的 曲率 连续 无 吸力 峰 前 缘 在 任意 工 况 
面 层 


a 


Ml 


总 可 以 消除 吸力 面 或 压力 面 近 前 缘 的 吸力 峰 ， 从 而 使 得 叶 表 阶 的 起 始 状 态 优 于 其 它 前 缘 叶 型 ， 因 此 气动 性 能 要 更 
好 。 

关键 词 : 叶 型 前 缘 ;， 亚 音 叶 型 ， 附 面 层 ， 数 值 分 析 ;， 压 气 机 
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Effects of Leading-Edge Geometry on Boundary Layer Parameters 


in Subsonic airfoil 


YU Xian-Jun, ZHU Hong-Wei, LIU Bao-Jie 
(National Key Laboratory of Science and Technology on Aero-Engine, Advanced Aero-Engine Collaboration Innovation Center, School 
of Energy and Power Engineering, Beihang University Beijing 100191, China) 
Abstract: The effect of circular, elliptical and curvature continuous leading edge on the development of boundary 
layer is studied by numerical Simulation method. It is found that the influence of the leading edge suction peak on the 
aerodynamic performance in subsonic compressor blade is due to the fact the suction affects the initial development 
state of the boundary layer by extracting the parameters of the boundary layer under different working conditions 
and the velocity pattern of the bounder layer. The leading edge of the suction peak diffusion parameter Dspike is too 
Strong will lead to the initial state of deterioration of the boundary layer, resulting in early transition, or even 
separation bubble, blade surface layer of rapid thickening, loss of blade increases, the available attack angle range 
reduced. The study also shows that the design of curvature continuous spikeless leading edge in any condition can 
always eliminate the Suction Surface or pressure surface near the leading edge of the suction peak, so that the initial 
state of the blade surface layer is better than the other leading edge, so aerodynamic performance should be better. 
Key words: Leading edge; Subsonic compressor blade; Boundary layer; Numerical analysis; Compressor 
0 引言 厚度 只 有 弱 长 的 1% 左 右 ， 在 工程 设计 中 往往 从 结构 
本 强度 层面 考虑 较 多 ; 另外 ， 在 叶 型 设计 过 程 中 ， 由 
未 来 飞机 推进 系统 要 求 压气 机 的 压 比 和 效率 都 ”于 前 尾 缘 曲 率 变化 较为 剧烈 ， 通 党 采用 分 段 设计 的 
大 大 高 于 目前 的 使 用 水 平 ， 以 提高 发 动机 的 推 重 比 ， ”方法 ， 即 先生 成 叶 身 再 添加 前 尾 缘 的 方式 ， 这 往往 
减 小 燃油 消耗 率 。 将 这 种 要 求 转换 到 叶片 上 即 为 : 会 导致 前 缘 和 叶 身 连接 处 的 曲率 变化 很 大 。 而 且 由 
高 负荷 、 大 可 用 攻 角 范围 和 健康 的 叶 表 附 面 层 趾 。 从 ” ”于 叶片 加 工 精度 不 足以 及 叶片 在 使 用 过 程 中 自然 侵 


国内 外 研究 趋势 上 看 ， 可 控 扩 散 叶 型 (CDA) 作为 。” 人 刨 和 磨损 等 原因 都 会 影响 叶片 前 缘 的 几何 形状 ， 这 
高 效率 、 高 负荷 压气 机 设计 的 关键 技术 之 一 ， 并 且 ” 些 几 何 形状 的 改变 都 会 对 叶片 性 能 产生 一 定 的 影响 。 
在 高 压 压气 机 中 后 部 的 亚 音 环境 中 己 经 得 到 了 非常 随 着 对 叶 型 的 深入 研究 ， 国 内 外 诸多 学 者 发 现 
广泛 的 工程 应 用 外 ,可 控 扩散 叶 型 设计 的 核心 思想 在 。”” 叶 型 前 缘 形 状 对 叶片 的 气动 性 能 有 着 不 可 忽视 的 影 
于 控制 叶 型 表面 的 等 炉 马赫 数 / 叶 表 静 压 升 系数 按照 ” 响 B， 随 之 这 也 成 为 了 叶 型 研究 的 热点 之 一 。 为 了 便 


特定 的 规律 分 布 ， 使 得 叶 表 附 面 层 发 展 的 更 加 健康 ， 于 加 工 ， 传 统 的 叶 型 前 缘 采 用 圆 形 设计 ， 但 是 随 着 
以 提高 叶 型 性 能 。 由 于 在 航空 发 动机 中 ， 压 气 机 叶 诸多 学 者 大 量 的 数值 模拟 和 实验 研究 表 发 现 *77， 由 
片 前 缘 厚 度 往往 很 小 ， 最 小 可 达 0.2mm 左右 ， 相 对 ”于 圆 形 前 缘 和 叶 身 连接 处 曲率 不 连续 导致 在 前 缘 附 


会 产生 分 离 泡 ， 造 成 流动 提前 转换 ， 增 加 损失 。 
陆 宏志 等 上 对 覃 圆 形 前 缘 进 行 了 数值 模拟 发 现 ， 椭 
圆 形 前 缘 能 够 明显 的 抑制 吸力 峰 ， 改 善 叶片 的 气动 
性 能 ， 作 者 基于 控制 前 缘 和 叶 身 曲率 过 度 的 想法 ， 
提出 了 带 平台 的 圆 弧 形 前 缘 ， 发 现 这 种 设计 可 以 改 
善 叶 表 流 动 ， 达 到 与 椭圆 形 前 缘 相 近 的 结果 。 鉴 于 
前 缘 问 题 的 存在 ， 宋 寅 &29 和 袁 春 香 00 等 人 又 提出 了 
前 缘 曲 率 连 续 的 设计 方法 ， 研 究 表明 ， 曲 率 连 续 的 
设计 方法 能 够 更 好 的 抑制 前 缘 吸 力 峰 〈Spike) 的 产 
生 ， 叶 表 附 面 层 的 发 展 更 加 健康 ， 很 有 效 的 改善 了 
叶片 性 能 。Andrew P. S Wheeler 和 MillerI12 等 对 上 
游 转子 尾 迹 叶 型 损失 的 影响 进行 的 实验 研究 表明 在 
上 游 尾 迹 干扰 的 情况 下 ， 椭 圆 形 前 缘 较 圆 形 前 缘 能 
够 明显 减 小 叶片 的 型 面 损失 。Goodhand 和 Miller03l 
对 椭圆 形 前 缘 进 行 优 化 消除 了 吸力 面前 缘 的 吸力 峰 
实现 了 设计 状态 无 前 缘 吸 力 峰 (spikeless ) 叶 型 设计 ; 
随 之 又 提出 吸力 峰 的 高 度 存 在 一 个 临界 值 ， 该 临界 
值 可 以 用 前 缘 吸 力 峰 扩散 因子 Peie 来 评价 ， 当 吸力 
峰 的 高 度 超 过 临界 值 后 ， 叶 型 损失 会 迅速 增加 。 
以 上 研究 结果 表明 ， 前 缘 形 状 对 叶 型 性 能 有 很 
大 的 影响 。 在 同一 状态 下 ， 不 同 的 前 缘 形 状 对 叶 型 
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Camber angle/(°) 41 
Turbulent intensity/(%) 3 
Reynolds number 


5.45*105 
(based on inlet velocity and chord) 
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1 前 缘 几 何 形状 
Figure 1 Leading-Edge geometry 
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性 能 的 影响 有 所 不 同 ， 但 是 前 缘 形 状 对 叶 表 附 面 层 
发 展 影响 的 机 制 还 没有 完全 澄清 。 鉴 于 此 ， 本 文 的 
主要 工作 就 是 针对 一 个 CDA 叶 型 ,在 不 改变 叶 身 的 
情况 下 分 别 实现 该 叶 型 的 圆 形 前 缘 、 椭 圆 形 前 缘 和 
Spikeless 前 缘 设 计 ， 并 利用 可 模拟 边界 层 转换 的 数 
值 分 析 方 法 对 叶 型 的 攻 角 损失 特性 进行 了 模拟 ， 然 
后 提取 了 详细 的 叶 表 附 面 层 在 的 发 展演 化 过 程 ， 对 
附 面 层 参 数 进 行 了 深入 分 析 。 通 过 对 比 不 同 前 缘 叶 
型 附 面 层 的 发 展 过 程 ， 探 讨 了 叶 型 前 缘 对 叶 型 气动 
性 能 影响 的 机 制 。 该 研究 将 为 进一步 的 实验 研究 和 
叶 型 优化 提供 必要 的 理论 支持 。 
1 数值 计算 方法 
1.1 叶 型 简介 

本 文采 用 的 可 控 扩 散 叶 型 参数 见 表 1。 叶 型 前 缘 
几何 形状 如 图 1。 


表 1 CDA 叶 型 参数 简介 
Table 1 CDA parameters 


Chord/m 0.18 
Solidity 1.6 
Inlet Mach number 0.13 
Stagger angle/(®) 21.77 
Inlet metal angle/(?) 47 


2 网 格 
Figure 2 mesh 
1.2 网 格 
本 文 主要 针对 附 面 层 流动 细节 进行 研究 ， 因 此 
对 网 格 要 求 较 高 。 网 格 采用 ICEM 生成 ， 在 叶片 表 
面 利 用 “OO” 型 网 格 在 保证 良好 的 正 交 性 的 同时 进行 
了 局 部 加 密 ， 保 证 了 近 壁 面 网 格 y+<1， 边 界 层 内 网 
ee 10 道 ， 且 网 格 以 1.1 的 比率 向 主流 区 递 
。 图 2 为 本 文采 用 的 网 格 ， 计 算 域 距 叶 片 前 缘 约 
了 倍 弦 长 ， 距 叶片 尾 缘 2 倍 弦 长 。 展 向 给 定 20 道 
eS 模型 网 格 总 数 为 64 万 。 
3 满 流 模型 选择 及 边界 条 件 设 定 
本 文 使 用 ANSYS CFX 软件 进行 叶 型 计算 ， 为 
了 更 加 准确 的 预测 附 面 层 的 转换 和 附 面 层 的 流动 细 
节 ， 计算 中 汕 流 模型 选取 Shear Stress Transport 模型 ， 
并 使 用 Gamma-theta 转换 模型 , 计算 格式 采用 高 精度 
格式 ， 保 证 二 阶 精 度 。 边 界 条 件 设 定 方式 为 : 进口 
边 条 给 定 总 温 、 总 压 和 气流 角 ; 出 口 边 条 根据 不 同 


的 工 况 给 定 出 口 流量 ; 栅 距 方向 设置 为 周期 性 边界 
条 件 ， 展 向 上 采用 对 称 边 条 ; 叶片 表面 为 光滑 无 滑 


移 固 体 表面 。 当 残 差 小 于 1.0X104 时 ， 判 定 计算 收 0.06 


敛 ， 但 需 保 证 流量 、 速 度 和 压力 残 差 曲线 平 直 。 一 Circle 


—o— Ellipse 
一 一 Spikeless 
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1006 0.2 0.4 0.6 0.8 1.0 Figure 4 Comparison of total pressure loss 
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3 叶 表 静 压 升 系数 45 


Figure 3 Coefficient of static pressure of blade Circle 


1.4 数值 模拟 校 验 


在 确定 了 数值 模拟 的 网 格 、 数 值 格式 、 江 流 模 注 ” 
型 、 转 换 模型 和 边界 条 件 之 后 ， 用 实验 数据 对 CFX 30 
的 计算 结果 进行 校 验 ， 以 确认 计算 结果 的 可 靠 性 。 
图 3 为 CFX 计算 结果 和 叶 栅 实验 结果 的 对 比 。 由 于 
在 叶 栅 实验 中 存在 定 的 三 维 流动 效应 ， 由 于 端 壁 200 35 34 36 38 40 42 44 46 48 50 52 54 
这 附 面 层 的 存在 ,进出 口 轴 流 密 流 比 AVDR 约 为 1.02， 和 
~。 ”但 是 在 数值 模拟 中 CFX 没有 考虑 AVDR 的 影响 ， 因 图 5 时 型 气流 折 转 角 比 较 
= 此 计算 结果 和 实验 结果 之 间 稍 有 偏差 ， 但 是 可 以 看 Figure 5 Comparison of flow deflection angle 
SI 到 计算 与 实验 结果 的 总 体 符合 性 较 好 ， 尤 其 是 在 前 


到 缘 附 近 的 计算 结果 和 实验 结果 符合 的 较 好 ， 对 吸力 
一 峰 的 模拟 与 实验 对 比 有 较 好 的 一 致 性 。 男 外 ， 由 于 
实验 中 没有 看 到 分 离 泡 型 附 面 层 转换 《可 能 是 由 于 
叶 型 的 叶 表 粗糙 度 以 及 来 流 消 流 度 稍 高 造成 的 )， 而 
个 计算 结果 显示 出 了 分 离 泡 转换 〈40% 弦 长 位 置 附近 ) 
〇 ”造成 了 压力 平台 。 考 虑 到 本 文 主要 分 析 的 是 前 缘 对 
叶 表 附 面 层 的 定性 有 影响， 因此， 数值 模拟 结果 具有 


0.2 
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较 好 的 参考 性 。 Inlet angle oy” 
a ey 图 6 叶 型 D 因子 比较 
2 本 缘 对 叶 型 气动 性 有 的 尿 Z 啊 Figure 6 Comparison of D factor 
图 4 给 出 了 计算 得 到 的 不 同 前 缘 几 何 造型 的 叶 rg 
型 总 压 损失 特性 ， 可 以 看 出 ， 三 种 前 缘 几 何 叶 型 的 Ellipse 


最 小 损失 基本 一 致 ， 但 是 可 用 攻 角 范围 《以 最 小 损 一 一 Spikeless 
失 的 两 倍 作 为 攻 角 计算 范围 ) 却 有 比较 大 的 差别 。 
圆 形 叶 型 的 可 用 攻 角 范围 明显 比 椭圆 和 Spikeless 前 
缘 叶 型 小 接近 2.5"， 接 近 总 工作 范围 的 14.3%。 
Spikeless 前 缘 时 型 可 用 攻 角 范围 最 宽 , 损失 也 最 小 。 


这 表明 不 同 的 前 缘 几 何 造成 了 叶 表 附 面 层 流 动 状态 0 32 3436 38 40 42 44 46 46 50 52 54 
的 差异 ， 导 致 了 损失 和 可 用 攻 角 范围 具有 差异 。 Inletiangle: em 


7 叶 型 落后 角 比 较 


Figure 7 Comparison of blade deviation angle 


从 图 5〈 叶 型 气流 折 转 角 特性 )、 图 6 ( 叶 型 D 
因子 特性 ) 和 图 7《〈 叶 型 落后 角 特 性 ) 中 可 以 看 到 ， 
圆 型 前 缘 叶 型 的 加 功 增 压 性 能 在 大 攻 角 状态 明显 要 
比 其 他 两 种 叶 型 差 , 椭圆 形 前 缘 叶 型 次 之 ,Spikeless 
前 缘 叶 型 的 性 能 最 优 。 还 可 以 看 到 ， 圆 形 前 缘 叶 型 
在 攻 角 范围 内 的 落后 角 都 明显 偏 高 ， 并 且 在 大 正 、 
负 攻 角 状 态 下 附 面 层 提前 分 离 ， 使 得 叶 型 损失 随 着 
攻 角 的 偏离 而 增加 。 而 Spikeless 前 缘 叶 型 落后 角 最 
小 ， 大 正 、 负 攻 角 状态 下 附 面 层 分 离 最 迟 ， 保 持 了 
很 好 的 流动 状态 。 综 上 所 述 ， 可 以 发 现 ， 在 其 他 参 
数 保持 一 致 ， 只 有 前 缘 几 何不 同 的 情况 下 ， 叶 型 性 
能 却 有 很 大 的 差别 ， 这 是 由 于 前 缘 几 何不 同 造 成 了 
叶 表 附 面 层 的 发 展 状态 出 现 差 异 ， 因 此 为 了 深入 而 
究 前 缘 几 何 对 叶 型 性 能 的 影响 ， 下 面 对 叶 表 附 面 层 
的 发 展 进行 详细 对 比分 析 。 


3 前 缘 对 叶 表 附 面 层 发 展 状态 的 影响 


3.1 附 面 层 参 数 提取 方法 

为 了 提取 附 面 层 参数 ， 自 主 发展 了 附 面 层 提取 
程序 ， 首 先 根据 计算 得 到 的 流 场 将 流 场 差 值 到 与 叶 
型 表面 严格 正 交 的 网 格 上 ， 然 后 在 正 交 网 格 上 进行 


附 面 层 参 数 的 提取 分 析 。 对 于 任意 位 置 的 附 面 层 厚 
度 利用 下 式 计算 : 


5=f (1 人 ee)dy GD 


pz2U 
积分 方向 沿 叶 表 法 向 , p14 为 当地 密度 和 速度 ; p20 为 
主流 密度 和 速度 。 其 中 主流 的 判断 根据 速度 梯度 沿 
法 向 的 极 值 点 来 判断 ， 当 没有 极 值 点 时 ， 则 根据 设 
定 梯 度 的 一 个 最 小 值 来 判断 。 由 于 边界 层 参 数 对 梯 
度 的 阀 值 不 敏感 ， 因 此 可 根据 主流 区 的 梯度 极 大 值 
适当 提高 作为 流 场 整体 的 边界 层 与 主流 的 分 界 判断 


准则 。 动 量 厚度 利用 下 式 计 算 : 
人 = 和 (1 一 各) dy (2) 
能 量 厚度 利用 下 式 计算 : 
"= (1 ei) dy (3) 
从 而 形状 因子 为 : 
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8 不 同 来 流 气 流 角 时 吸力 面 附 面 层 形状 因子 的 发 展 过 程 ， 
来 流 气 流 角 (a)35°(b)47° (c)50° 
Figure 8 Development of boundary layer shape factor on suction 


surface at inflow angle of (a)35°(b)47° (c)50° 


3.2 不 同 攻 角 状态 叶 型 附 面 层 发 展 规 律 分 析 

图 8、 图 9 和 图 10 为 不 同 攻 角 下 叶 型 吸力 面 附 
面 层 积分 参数 的 变化 ， 从 中 可 以 看 出 ， 在 来 流 气 流 
角 为 35° 的 大 负 攻 角 工 况 下 ,三 种 叶 型 的 滞 止 点 位 置 
均 位 于 吸力 面 一 侧 ， 在 吸力 面 一 侧 的 前 缘 绕 流 过 程 
没有 很 大 的 差别 ， 因 此 附 面 层 发 展 过 程 基本 一 致 ， 
都 是 在 60% 流 向 位 置 处 发 生 转换 ， 导 致 附 面 层 迅速 
增 厚 ， 摩 阻 损 失 增 大 ， 因 此 在 吸力 面 的 损失 基本 
致 。 在 来 流 气流 角 为 47° 的 设计 工 况 下 ， 三 种 叶 型 的 
灌 止 点 位 置 均 在 前 缘 点 附近 ， 由 于 Spikeless 叶 型 前 
缘 曲 率 变 化 最 小 ，Ellipse 叶 型 前 缘 曲 率 变化 次 之 ， 
两 者 在 前 缘 吸 力 面 一 侧 的 绕 流 过 程 很 相近 ， 都 没有 
发 生 过 强 的 突然 压缩 和 膨胀 过 程 ， 都 是 在 40% 流 向 
位 置 发 生 了 转换 ， 造 成 附 面 层 迅速 增 厚 ， 摩 阻 损失 
增 大 。 但 是 circle 叶 型 前 缘 处 存在 曲率 不 连续 ， 从 而 
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circle 叶 型 的 前 缘 绕 流 存在 较 强 的 突然 压缩 和 过 度 脱 ee 

胀 过 程 , 即 较 强 的 前 缘 Spike, 并 在 前 缘 出 现 分 离 泡 ， Ellipse 
导致 吸力 面 在 前 缘 附 近 就 发 生 分 离 泡 转换 ， 附 面 层 — 人 — Spikeless 
也 随 之 迅速 增 厚 ， 摩 阻 损失 增 大 。 因 此 ， 在 该 状态 
下 circle 叶 型 吸力 面 的 损失 最 大 , Ellipse 和 Spikeless 
的 损失 基本 一 致 且 最 小 。 
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10 不 同 来 流 气 流 角 时 吸力 面 附 面 层 表面 摩 阻 系数 的 发 展 ， 
来 流 气 流 角 (a)35°(b)47°(c)50° 
04 06 08 10 Figure 10 Development of boundary layer friction coefficient 


| 
on Suction surface at inflow angle of (a)35°(b)47° (c)50° 

在 来 流 气 流 角 为 50° 的 大 正 攻 角 工 况 下 , 三 种 叶 
型 的 灌 止 点 位 置 均 位 于 压力 面 一 侧 ， 流 体 需要 绕 过 
前 缘 点 到 达 吸 力 面 ， 而 三 种 叶 型 前 缘 的 曲率 变化 有 
较 大 的 差别 ， 这 将 导致 前 缘 绕 流 存在 不 同 程度 压缩 
和 膛 胀 过 程 。 由 于 circle 叶 型 前 缘 除 前 缘 点 外 曲率 最 
大 并 保持 一 致 ， 但 是 在 前 缘 和 叶 身 连接 点 存在 曲率 


00 02 04 05 08 10 突变 , 因此 circle 叶 型 的 前 缘 绕 流 的 突然 压缩 和 过 度 
9 Cc) 膨胀 过 程 最 强 ， 前 缘 Spike 也 最 强 , 并 在 吸力 面前 缘 

9 不 同 来 流 气 流 角 时 吸力 面 附 面 层 能 量 厚度 的 发 展 过 程 ， 出 现 分 离 泡 。 吸 力 面 转换 位 置 提前 到 了 前 缘 附近 ， 
来 流 气流 角 ()35°(b)47*(c)50* 造成 附 面 层 过 快 的 增 厚 ， 并 在 尾 缘 出 现 了 较 大 的 附 
Figure 9 Development of boundary layer energy thickness on 面 层 分 离 ， 增 加 了 吸力 面 叶 型 损失 。 Ellipse 叶 型 前 
suction surface at inflow angle of (035(b)47 (Oj50” 。。。 缘 点 曲率 较 大 ， 除 前 缘 点 外 前 缘 其 他 位 置 曲率 较 小 ， 
一 = 一 Circle 而 且 前 缘 和 叶 身 连接 处 的 曲率 突变 较 小 ,因此 Ellipse 
| 一 一 Ellipse 叶 型 前 缘 绕 流 造成 的 突然 压缩 和 过 度 膨 胀 过 程 较 


0020 上 全 - Spikeless Circle 前 缘 弱 ， 虽 然 转换 位 置 也 提前 到 了 前 缘 附 近 ， 
但 是 造成 的 附 面 层 增 厚 却 较 慢 ， 吸 力 面 叶 型 损失 增 
加 较 circle 叶 型 小 。 虽 然 Spikeless 叶 型 前 缘 点 曲率 
最 大 ， 但 是 除 前 缘 点 外 前 缘 其 他 位 置 曲 率 最 小 ， 前 
s | 缘 和 叶 身 曲率 连续 ， 这 大 大 减弱 了 前 缘 绕 流 的 膨胀 
0 0 06 “08 10o 过 程 ， 优 化 了 前 缘 绕 流 过 程 ， 叶 型 吸力 面 转换 位 置 

X/Xo (a) 依然 保持 在 30% 的 流向 位 置 , 比 其 他 两 种 叶 型 靠 后 ， 
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而 且 尾 缘 的 附 面 层 分 离 也 最 弱 ， 损 失 也 最 小 。 
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11 不 同 来 流 气流 角 时 压力 面 附 面 层 形状 因子 的 发 展 过 程 ， 
来 流 气流 角 (a)35°(b)47°(c)50° 
Figure 11 Development of boundary layer shape factor on 
pressure surface at inflow angle of (a)35°(b)47° (c)50° 
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12 不 同 来 流 气流 角 时 压力 面 附 面 层 能 量 厚度 的 发 展 过 程 ， 
来 流 气 流 角 (a)35°(b)47°(c)50° 
Figure 12 Development of boundary layer energy thickness on 


pressure surface at inflow angle of (a)35°(b)47° (c)50° 


图 11、 图 12 和 图 13 为 不 同 攻 角 下 叶 型 压力 面 
附 面 层 积分 参数 的 变化 。 可 以 看 到 ， 在 来 流 气 流 角 
为 3" 的 大 负 攻 角 工 况 下 , 三 种 叶 型 的 潇 止 点 位 于 吸 
力 面 一 侧 ， 因 此 前 缘 绕 流 过 程 和 上 文 所 述 的 大 正 攻 
角 工 况 下 吸力 面 的 绕 流 过 程 类 似 ， 虽 然 压力 面 转换 
都 在 前 缘 发 生 〈 分 离 泡 转换 )， 但 是 circle 叶 型 的 前 
缘 绕 流 的 突然 压缩 和 膨胀 过 程 最 强 , ellipse 叶 型 次 之 ， 
Spikeless 最 弱 。 这 就 造成 Circle 压力 面 附 面 层 增 厚 最 
快 ， 厚 度 最 大 ，Ellipse 叶 型 次 之 ，Spikeless 叶 型 最 
小 。 因 此 , 在 该 状态 下 , Circle 叶 型 压力 面 损失 最 大 ， 
Ellipse 叶 型 其 次 ，Spikeless 叶 型 最 小 。 
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13 不 同 来 流 气流 角 时 吸力 面 附 面 层 表面 摩 阻 系数 的 发 展 ， 
来 流 气流 角 (a)35°(b)47°(c)50° 
Figure 13 Development of boundary layer friction coefficient on 
pressure surface at inflow angle of (a)35°(b)47° (c)50° 

在 来 流 气流 角 为 47" 的 设计 工 况 下 ,三 种 叶 型 的 
滞 目 点 均 位 于 前 缘 点 附近 ， 由 于 在 吸力 面 绕 流 的 压 
缩 和 膨胀 过 程 较 强 ， 因 此 在 压力 面 的 压缩 和 膨胀 过 
程 相应 就 会 弱 ， 三 种 叶 型 均 没 有 发 生前 缘 压 力 面 转 
换 《〈 形 状 因子 在 2.5 左右 )， 此 时 压力 面 没 有 发 生 分 
离 泡 转换 。 从 形状 因子 来 看 ， 一 直到 尾 缘 ， 压 力 面 
的 形状 因子 都 保持 在 2 以 上 ， 因 此 可 以 判断 此 时 模 
拟 得 到 的 压力 面 附 面 层 状态 属于 全 层 流 状态 ， 故 总 
体 上 损失 都 较 小 。 但 是 不 难看 出 ,Ellipse 前 缘 和 Circle 
前 缘 叶 型 由 于 前 缘 与 叶 身 之 间 曲 率 过 度 存在 一 定 的 
不 连续 性 ， 使 得 在 靠近 前 缘 区 域 存在 一 个 局 部 扰动 ， 
但 没有 触发 典型 的 层 流 边 界 层 分 析 ， 并 且 在 该 区 域 
摩 阻 系数 存在 一 个 尖峰 。 由 于 前 缘 区 域 对 附 面 层 造 
成 的 抠 动 不 同 ， 因 此 附 面 层 的 发 展 起 始 状态 也 有 所 
不 同 ， 导 致 了 从 前 缘 到 尾 缘 附 面 层 的 发 展 状态 还 是 
不 同 ， 但 是 三 者 损失 相差 不 大 。 当 来 流 气 流 角 增 大 
到 50° 的 大 攻 角 工 况 时 ,压力 面 一 侧 前 缘 绕 流 的 压缩 
和 膨胀 过 程 进 一 步 减弱 ， 三 种 叶 型 的 压力 面 都 保持 
为 层 流 状 态 ， 流 动情 况 与 47° 工 况 相 似 。 
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图 14 吸力 面 不 同 流向 位 置 附 面 层 速度 型 〈 来 流 气 流 角 50°) 
Figure 14 Boundary layer velocity profile on blade suction 
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surface at inflow angle of 50° 


由 以 上 分 析 可 以 看 到 ， 前 缘 不 同 带 来 的 前 缘 绕 
流 过 程 会 有 所 差异 ， 从 而 引起 附 面 层 的 厚度 不 同 ， 


损失 也 会 不 同 。 如 图 14 所 示 ( 模 坐标 V 表示 速度 大 
小 ， 纵 坐标 sh 表示 距离 当地 叶片 表面 的 法 向 高 度 )， 
通过 观察 不 同 的 流向 位 置 附 面 层 速度 型 的 分 布 ， 可 
以 更 为 清晰 地 分 析 附 面 层 速度 剖面 的 发 展 过 程 。 从 
图 中 可 以 看 到 ， 在 来 流 角度 为 50° 时 ，Circle 叶 型 前 


缘 绕 流 的 加 速 和 膨胀 过 程 最 严重 ， 从 X=0.0005 到 
X=0.0024 Circle 前 缘 时 型 的 边界 层 速度 剖面 从 最 饱 
满 变 成 了 近 壁 亏损 最 严重 ， 并 且 可 以 明显 看 到 此 时 
前 缘 的 分 离 泡 从 X=0.0009 前 已 经 起 始 , 在 X=0.0024 
之 后 结束 。 对 于 该 分 离 泡 的 存在 性 从 图 10 给 出 的 叶 
表 摩 阻 系数 分 布 也 可 以 看 到 ， 即 在 近 前 缘 区 域 存在 
一 个 Ct 突 降 为 零 的 波动 。 从 X=0.0075 的 结果 中 可 以 
看 出 ，Circle 前 缘 叶 型 的 前 缘分 离 泡 已 经 再 附 , 但 是 
速度 型 的 饱满 程度 显然 较 男 外 两 个 叶 型 都 要 差 很 多 。 
图 8 形状 因子 的 分 布 上 来 看 ， 经 过 分 离 泡 之 后 附 
层 应 当 已 经 发 生 了 转换 ， 因 此 附 面 层 的 厚度 增加 
很 快 ， 明 显 高 于 另外 两 种 前 缘 叶 型 的 附 面 层 。 
从 图 14 所 示 的 Spikeless 前 缘 叶 型 的 吸力 面 附 面 
层 剖 面 的 发 展 过 程 可 以 看 出 ， 在 近 前 缘 区 域 没 有 分 
离 的 特征 ， 因 此 ， 没 有 出 现 分 离 泡 ， 这 与 图 10 所 示 
的 结果 表现 一 致 。 从 X=0.0009 之 后 ，Spikeless 前 缘 
叶 型 吸力 面 的 速度 剖面 在 三 个 叶 型 中 都 最 为 饱满 ， 
附 面 层 厚 度 明显 最 小 ,从 图 10 所 示 的 结果 可 以 看 出 ， 
来 流 角 度 为 50° 时 Spikeless 前 缘 叶 型 吸力 面 在 20~30% 
弦 长 之 间 的 摩 阻 系数 非常 接近 0, 此 时 附 面 层 应 当 是 
非常 接近 分 离 的 状态 ， 并 且 在 此 过 程 之 中 发 生 了 附 
面 层 的 转换 《〈 见 图 8 所 示 的 结果 )， 对 应 的 X 位 置 应 
当 在 X=0.036 到 0.054 之 间 。 从 图 14 所 示 的 结果 可 
以 看 到 ， 在 X=0.022 位 置 Spikeless 前 缘 叶 型 的 附 面 
层 的 近 壁 速度 梯度 明显 小 于 另外 两 种 叶 型 ， 表 现 为 
典型 的 层 流 速度 剖面 ， 但 是 在 X=0.06 的 位 置 已 经 转 
换 成 了 典型 的 油 流 速度 型 ， 说 明 在 这 个 过 程 中 发 生 
了 转换 ， 并 且 应 当 是 尺度 相对 较 小 的 分 离 泡 转换 。 
由 于 Spikeless 前 缘 时 型 转换 的 位 置 最 为 靠 后 ， 因 此 
后 续 发 展 的 满 流 附 面 层 也 相对 更 为 健康 。 

从 图 14 所 示 的 Ellipse 前 缘 叶 型 叶 表 附 面 层 剖 面 
的 发 展 情况 ， 可 以 看 出 Ellipse 前 缘 叶 型 的 前 缘分 离 
泡 尺 度 与 Circle 前 缘 叶 型 相 比 尺度 相当 ,但 是 强度 稍 
弱 ， 分 离 泡 之 后 也 发 生 了 完全 转换 ， 因 此 吸力 面 叶 
附 面 层 是 层 流 附 面 层 。 相 关 现 象 与 图 8 和 图 10 中 
显示 的 结果 一 致 。 总 体 上 ，Ellipse 前 缘 叶 型 的 附 面 
层 厚 度 和 饱满 性 介 于 Circle 前 缘 叶 型 和 Spikeless 前 
缘 叶 型 之 间 。 

综 上 所 述 ， 可 以 发 现 ， 即 使 在 前 缘 前 后 经 历 相 
同 的 逆 压 梯度 ， 不 同 的 前 缘 几 何 造成 的 前 缘 绕 流 过 
程 不 同 ， 直 接 导 致 附 面 层 的 起 始 发 展 状态 不 同 ， 使 
得 叶 表 流动 状态 产生 了 很 大 差异 ， 造 成 能 量 亏 损 不 
同 〈 如 图 9 和 图 12 所 示 )， 即 叶 型 损失 不 同 。 这 也 
就 造成 前 缘 几 何不 同 ， 叶 型 损失 有 很 大 的 差别 ， 可 
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用 攻 角 范围 也 会 有 很 大 的 变化 。 由 此 可 见 ， 不 同 的 
前 缘 形 状 对 叶片 性 能 的 影响 与 前 缘 绕 流 强度 息 息 相 
关 ， 该 绕 流 过 程 中 造成 的 局 部 加 速 和 减速 越 强 ， 通 
常 流动 损失 也 越 大 。 

4 前 缘 局 部 绕 流 对 叶 表 附 面 层 的 影响 
4.1 不 同 前 缘 造 型 对 前 缘 Spike 发 展 的 影响 

为 了 分 析 前 缘 局 部 扰 流 的 关键 影响 因素 ， 有 必 


要 借鉴 Goodhand 和 Miller03] 提 出 的 Dspixe 参数 进行 
深入 分 析 。 该 参数 的 定义 如 下 : 


_ UMAX—-UMIN 
Dspire = or (5) 


其 中 各 变量 的 定义 如 图 15 所 示 。 该 参数 实际 上 
是 借鉴 D 因子 提出 来 的 ， 表 征 了 前 缘 突 然 压缩 过 程 
的 强 弱 。 


Suction-surface spike S 


15 叶片 表面 压 升 系数 和 前 缘 局 部 压 升 系数 分 布 示意 图 53 
Figure 15 The static pressure coefficient of blade surface and the 
local static pressure coefficient of the leading edgel!3] 

图 16 为 不 同 前 缘 叶 型 叶 表 压 升 系数 随 来 流 气 流 
角 的 变化 ， 结 合 图 17 给 出 的 吸力 面前 缘 Dspixe 随 来 

流 气流 角 的 变化 , 可 以 看 出 , 随 着 来 流 气流 角 增 大 ， 
吸力 面前 缘 Dspixe 也 随 之 增 大 ， 即 前 缘 吸 力 峰 的 强度 
增 大 。 此 外 ， 不 同 的 前 缘 几 何 叶 型 Dspie 增 长 速率 并 
不 一 致 ，Circle 前 缘 叶 型 Dsoike 增长 最 快 ， 更 容易 出 
现 吸力 面 分 离 泡 ， 增 加 吸力 面 叶 型 损失 。Spikeless 
叶 型 在 负 攻 角 状 态 保持 很 低 的 Dspxe， 并 且 随 来 流 气 
流 角 的 变化 增长 缓慢 ， 只 有 在 45? 攻 角 之 后 ，PDaspie 
增长 速率 加 快 。Ellipse 前 缘 叶 型 Dspixre 的 变化 介 于 
Circle 前 缘 叶 型 和 Spikeless 叶 型 之 间 。 由 图 18 看 出 ， 
相对 于 吸力 面 而 言 ， 压 力 面前 缘 绕 流 造 成 前 缘 Spike 
较 小 ， 随 着 来 流 气 流 角 的 增 大 ， 压 力 面 前 缘 Dspre 减 
小 ， 对 叶 型 性 能 的 影响 也 随 之 减 小 。 此 外 ， 可 以 明 
显 看 到 Spikeless 前 缘 的 设计 实现 了 吸力 面 从 叶 型 最 
小 损失 攻 角 状态 《〈45" 左 右 ) 到 所 有 负 攻 角 状 态 的 无 
吸力 峰 流 动 ， 而 压力 面 从 最 小 损失 攻 角 状态 到 正 攻 
角 边 界 之 间 无 吸力 峰 ， 即 所 设计 的 Spikeless 前 缘 总 
可 以 实现 吸力 面 或 者 压力 面 的 无 Spike 流动 。Circle 
前 缘 叶 型 除了 在 最 大 负 攻 角 状 态 吸 力 面 Spike 较 弱 ， 


最 大 正 攻 角 状态 压力 面 Spike 较 弱 外 ， 其 它 状 态 的 
Spike 都 非常 明显 。Ellipse 前 缘 在 负 攻 角 边 界 附近 实 
现 了 吸力 面 无 Spike, 最 大 正 攻 角 状态 压力 面 无 Spike， 
其 它 状态 都 存在 Spike。 显 然 ，Spike 的 强 弱 与 叶 表 
附 面 层 的 发 展 状态 直接 关联 ， 即 与 叶 型 损失 直接 相 
关 。Spikeless 前 缘 叶 型 之 所 以 能 够 实现 好 的 叶 型 性 
能 , 与 基本 消除 或 者 显著 削弱 了 Spike 密切 关联 ， 并 


1 . 1 2. 1 . i 2 
35 40 45 50 55 


且 Spikeless 的 状态 较 宽广 。 因 此 ， 有 必要 对 叶 表 附 Inlet angle of/° 
面 层 参 数 与 Spike 之 间 的 关联 进行 深入 分 析 。 图 17 吸力 面前 缘 Dwae 随 来 流 气 流 角 的 变化 
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Figure 17 Variations of Dspike with the inlet angle at the 
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4.2 前 缘 Spike 与 叶 表 附 面 层 参 数 之 间 的 关系 


一 时 一 Inlet angle 35° 
一 C 一 Inlet angle 47° 
一 Inlet angle 50° 


0.00 005 ~ 010 0.15 0.20 基于 上 面 的 分 析 ， 提 取 了 10% 弦 长 位 置 的 叶 表 
X/X, (b) 


动量 厚度 和 形状 因子 ,分 析 Dspike 与 附 面 层 参 数 之 间 
SR 的 关系 。 从 图 19 和 图 20 可 以 看 出 ， 叶 型 能 量 损失 
一 2 一 Inletangle47* 厚度 随 着 Dspixe 的 增 大 而 增 大 ， 当 Desiee 达到 一 定 值 
之 后 由 于 分 离 泡 的 出 现 ， 导 致 了 叶 型 损失 迅速 增 大 。 
对 于 图 19 所 示 的 吸力 面 结 果 ，Spikeless 前 缘 叶 型 与 
其 它 前 缘 叶 型 相 比 在 Dspire 相同 的 情况 下 ， 造 成 的 能 
量 损失 厚度 要 更 大 一 些 ，Ellipse 前 缘 叶 型 与 Circle 
前 缘 叶 型 相 比 也 要 大 一 些 。 此 外 ， 从 图 19 所 示 的 结 
可 以 看 到 Dspire>0.17 后 才 会 造成 吸力 面 叶 型 损 
失 的 显著 增 大 (Spikeless 前 缘 和 Circle 前 缘 叶 型 的 临 
界 值 接近 0.2)， 这 与 Goodhand 和 Miller03 给 出 的 
Dasoie>0.1 的 准则 有 较 大 偏差 。 而 从 图 20 所 示 的 压力 
结果 可 以 看 出 ， 此 时 不 同 前 缘 叶 型 Dspixe 与 附 面 层 
能 量 损失 厚度 的 关系 非常 接近 ， 规 律 与 吸力 面 也 相 
同 ，Circle 前 缘 叶 型 在 相同 的 Dspixe 情况 下 造成 的 损 
失 要 稍 小 一 些 。 此 外 ， 压 力 面 的 结果 显示 Dspire>0.12 
时 会 造成 压力 面 叶 型 损失 的 显著 增 大 (Circle 前 缘 叶 
型 的 临界 值 接近 0.14)， 这 与 Goodhand 和 Miller3l 
给 出 的 准则 较为 接近 ; 但 是 需要 指出 的 是 ， 此 时 
Circle 前 缘 叶 型 叶 型 的 临界 值 接 近 0.18。 由 此 见 ， 从 
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Figure 16 Variations of the local static pressure coefficient near 
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Figure 19 Variation of boundary layer energy thickness with 


Dspike on suction surface at 10% chord location 
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Figure 20 Variation of boundary layer energy thickness with 


Dspike on pressure surface at 10% chord location 
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Figure 21 Variation of boundary layer shape factor with Dspike 


on suction surface at 10% chord location 
图 21 和 图 22 中 给 出 了 叶 型 前 缘 局 部 区 域 吸力 
和 压力 面 边界 层 形状 因子 与 Dspire 的 关联 关系 。 可 
以 看 出 ,在 吸力 面 Dspixe>0.17 后 附 面 层 的 形状 因子 迅 
速 减 小 ， 从 层 流 状态 转变 为 清流 状态 ， 即 发 生 了 转 
换 。 这 表明 通常 情况 下 吸力 面 Dspire>0.17 会 造成 前 缘 
区 域 出 现 分 离 泡 并 触发 转换 。 同 样 可 以 看 到 ， 同 在 


SS 


层 流 状态 ，Circle 前 缘 叶 型 的 形状 因子 要 小 一 些 ， 即 
附 面 层 剖 面 的 近 壁 速度 梯度 要 更 小 一 些 ， 这 是 造成 
能 量 损失 厚度 要 小 一 些 的 原因 ， 这 应 当 与 前 缘 局 部 
的 曲率 有 关 。 从 图 22 显示 的 压力 面 结 果 可 以 看 出 ， 
此 时 诱发 边界 层 转换 的 临界 Dspie 约 为 0.12〈Circle 
前 缘 叶 型 的 临界 值 接近 0.14)， 由 于 形状 因子 迅速 减 
小 后 有 快速 增 大 ， 这 说 明 如 果 Dspire 如 果 过 强 
(Spikeless 前 缘 和 Ellipse 前 缘 叶 型 Dspike>0.13,Circle 
前 缘 叶 型 Dspixe>0.15) 极 有 可 能 直接 造成 了 大 尺度 的 
负 攻 角 压 力 面 分 离 流动 ， 损 失 会 快速 增加 。 
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Figure 22 Variation of boundary layer shape factor with Dspike on 


pressure surface at 10% chord location 
5 结论 


本 文 以 三 种 不 同 前 缘 造 型 的 叶 型 为 例 ， 通 过 数 
值 模拟 研究 了 不 同 前 缘 造 型 对 叶 型 表面 附 面 层 发 展 
机 制 ， 得 到 以 下 主要 如 下 结论 : 

(1)、 不 同 的 前 缘 造 型 会 使 得 前 缘 绕 流 过 程 不 
同 ， 即 前 缘 Spike 不 同 ， 过 强 的 前 缘 Spike 会 导致 叶 

4 性 能 下 降 ， 可 用 攻 角 范围 减 小 。Circle 叶 型 前 缘 绕 
流 过 程 变 化 最 为 剧烈 ， 损 失 偏 大 ， 可 用 攻 角 范围 较 
小 ，Spikeless 叶 型 前 弱 了 前 缘 Spike， 减 小 了 叶 型 损 
失 ， 增 大 了 可 用 攻 角 范围 。 

(2)、 前 缘 Spike 对 叶 型 性 能 影响 的 主要 在 于 ， 
前 缘 Spike 影响 了 附 面 层 的 起 始 发 展 状 态 , 过 强 的 前 
缘 Spike 使 得 速度 型 不 饱满 , 在 经 历 相同 逆 压 力 梯 度 
的 情况 下 ， 更 难以 抵抗 逆 压 梯度 ， 造 成 附 面 层 提前 
转换 ,甚至 出 现 分 离 泡 , 从 而 导致 附 面 层 迅速 增 厚 ， 
叶 型 损失 增加 ， 可 用 攻 角 范围 减 小 。 

(3)、 吸 力 面前 缘 Dspire 随 着 来 流 气 流 角 的 增 大 
而 增 大 ， 对 吸力 面 性 能 的 影响 也 随 着 增 大 ， 使 得 吸 
力 面 叶 型 性 能 下 降 ， 压 力 面 反之 。 

(4)、 叶 型 损失 随 着 前 缘 Dspire 的 增 大 而 增 大 ， 
当 达 到 某 一 值 后 ， 会 引起 叶 型 性 能 的 急剧 下 降 。 


际 


(5)、 所 设计 的 Spikeless 前 缘 在 任意 工 况 总 可 
以 消除 吸力 面 或 者 压力 面 近 前 缘 的 吸力 峰 ， 从 而 使 
得 叶 表 附 面 层 的 起 始 状 态 优 于 其 它 前 缘 叶 型 。 
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